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Abstract 



We deduce the gluino contribution to the three-loop QCD (3 function within 
the minimal supersymmetric Standard Model (MSSM) from its standard QCD 
expression. The result is a first step in the computation of the full MSSM 
three-loop (3 function. In addition, in the case of a light gluino it provides the 
strong three-loop SUSY correction to the extrapolation of the strong coupling 
constant from the low energy regime to the Z region and up to the squark 
threshold. 
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Typeset using REVTeX 



Although experimental measurements at the highest available energy are consistent with 
the standard model Q, the observed relationship of the strong coupling constant at the Z 
and the weak angle as well as the value of the b/r mass ratio vis-a-vis the top quark mass 
remain strong indications of a supersymmetric (SUSY) grand unification above 10 w GeV and 
a SUSY threshold for squarks and sleptons in the 0.1 to 1 TeV region. In this unification 
picture the value of the SUSY threshold is very sensitive to the highest known (two-loop) 
contribution to the MSSM (3 function. At the one-loop order a SUSY threshold far below 100 
GeV would be needed to fit the coupling constant measurements and such a low threshold 
is directly ruled out by the non-observation of squarks and sleptons in Z decay. This can be 
seen by inserting the current experimental values for the couplings into the semi-analytic ex- 
pressions for the SUSY scale given in . This suggests that the three-loop results could also 
be important especially as the precision of the measurements at the Z and beyond improves. 
As a first step in the calculation of the full three-loop (3 functions of the MSSM, we provide 
here the gluino contribution to the renormalization of the strong coupling constant. This 
gives the complete result in the region between the gluino mass and the squark mass which, 
in the light gluino scenario, extends from the low energy regime up to the Z and beyond up 
to the squark threshold. The standard Lagrangian density of the MSSM is summarized in 
||. The gluons interact with quarks and squarks in the fundamental representation of the 
gauge group SU(N), and with gluons, gluinos, and ghosts in the adjoint representation. In 
each representation the generators satisfy the commutation relations 

[R a ,R b ]=if abc R c , (1) 

with the adjoint representation matrices being defined in terms of the structure constants 
{F a )b c = if bac . The running of the strong coupling constant as a function of the scale /i is 
determined by the QCD (3 function 

da s /d(ln jj 2 ) = a s (3(a s / An), (2) 

where (3 has the perturbative expansion 

p(x) = Pix + (3 2 x 2 + p 3 x 3 + .... (3) 

Ignoring squark contributions, the one- and two-loop results in the minimally extended 
SUSY QCD B are 

A = ~^ + |(n/r+^), (4) 

A = ~<? A + f(n f TC A + f Cl) +A(n f TC F + ^C A ), (5) 

where n/ is the number of quark flavors and n g is the number of gluino multiplets, C A and 
Cf are the eigenvalues of the quadratic Casimir operators in the adjoint and fundamental 
representation respectively and T is the Dynkin index for the fundamental representation. 
Note that, throughout this paper we use the framework of the dimensional regularization || 
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and a minimal subtraction (MS, MS) type prescriptions ||. In the Standard Model (n g = 0), 
the three-loop coefficient is pfl 

0T = -^cj - »,r ( 2C| - ^c F c, - If c3) - („,iy ( «c P + ^). (e) 

Arriving at this result required the calculation of several hundred Feynman graphs, and a 
similar number is required to extend the result to the MSSM. However the first step in this 
program, the incorporation of gluino loops, can be obtained by purely group theoretical 
methods. This leaves a sharply reduced number of graphs that must be treated in detail 
to determine the full three loop f3 function. The main result of this paper, proven in the 
appendix, is that, in the MSSM above the gluino mass scale but below that of the squarks, 
there are only six linearly independent group factors, Cj, i — 0, ... ,5, among all the relevant 
Feynman graphs. In terms of these the three-loop (3 function of the MSSM takes the form: 

& = b C + hd + b 2 C 2 + b 3 C 3 + 64C4 + &5C5. (7) 

The graphs involving solely gluons and ghosts have each a group factor proportional to 
Co = C\ as can be trivially deduced on dimensional grounds from the case where there are 
no quarks in the theory. Graphs involving a single fermion loop have group weights that are 
linear combinations of the three factors 

C l = n f TC 2 F + n g C A /2, (8) 
C 2 = n f TC F C A + n g C 3 A /2, (9) 

C 3 = n f TC 2 A + n 9 C A /2. (10) 

Finally, all graphs involving two fermion loops have group weights that are linear combina- 
tions of the two factors 

C 4 = (n f TC F + n g C\/2){n f T + n g C A /2), (11) 

C 5 = (n f TC A + n 9 C 2 A /2)(n f T + n~ g C A /2). (12) 

In the standard model case where n g = 0, these results are trivial and not useful. The 
importance of eq. (§) - ( |T2"D is that they constrain the ways contributions from gluino loops 
follow from those of quark loops. Substituting Co through C5 into eq. (0) and comparing 
with eq. (||) in the n g = limit suffices to determine the coefficients bo through 65. The final 
result for the MSSM including gluinos but excluding squark contributions is therefore 

Ps - -^Ci - n f T { 2C F - 2 f CF C A - l -fcl) - { n f Tf^C F + f C A ) 

+ 9 -§n~ g C A - n g n f T( 2 -^C A + fc A C F ) - ^np A . (13) 
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Assuming the gluino lies below the Z and the squarks above the Z, then at the Z scale, f3$ 
is -9769/54 in the standard model (n s = 0) and +14134/27 in the MSSM (n s = 1). If the 
gluino lies below the b quark, the value of a s at the Z for a given a s at is increased by 
a non-negligible amount compared to the experimental error. A complete analysis of the 
precise effect is left for a later complete phenomenological analysis which should include 
the light gluino effect on the Z and r decay widths ||||. The results found here for the 
MSSM using the dimensional regularization || framework, as expected, do not match to 
those found recently [|iO| , |TT| within the dimensional reduction framework and they do not 
have to. A similar scheme dependence has been found in |]12|| . 
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APPENDIX: 

The group factors occuring in the f3 function are calculated in terms of traces of repre- 
sentation matrices. The fundamental representation matrices, T a , and the adjoint represen- 
tation matrices, F a satisfy 

T a T a = C F 1, Tr{T a T h ) = T5 ab , (Al) 

F a F a = C A 1, Tr{F a F b ) = C A 8 ab , (A2) 

papbpa = C A F b /2, T a T b T a = ^ _ C A /2)T b , (A3) 

Tr(T a T b T c ) = T{d abc + if abc ) /2, Tr(F a F b F c ) = if abc C A /2, (A4) 

j abc j abc = Ca ( N 2 _ ^ A5 ^ 

where 

C A = 2TN, C F = T(N 2 - 1)/N. (A6) 

The unit matrices appearing in eqs. (|A1|) , (|A2|) are N and N 2 — 1 dimensional for the fun- 
damental and adjoint representations of SU(N) respectively. The arbitrary normalization 
of the generators is usually chosen so that T = 1/2. The gauge coupling constant usually 
quoted in experimental analyses uses this normalization. The result of this work, that the 
three-loop (3 function including gluino but not squark contributions is a linear combination 
of the six quoted group factors, is proven in this appendix. 

We begin by noting that the MS scheme strong coupling renormalization constant Z as 
can be determined via the renormalization constants Z\ for the ghost-ghost-gluon vertex, 



4 



Z\j 2 for the gluon propagator, and Z l J 2 for the ghost propagator (see, e.g., [T^] for more 
details) 

Z as = ZlZ^Zf, (A7) 

where by the definition of the MS type schemes each renormalization constant Zi is a 
polynomial 



= l + (a s )e- n 7 (A8) 

n>l 



with e = (4 — D)/2. Now, the QCD (3 function can be obtained from the following relation 
(see, e.g., 0) 

dzw 

P = ^r^- (A9) 

am a s 

Contributions from graphs with no fermion loops each have the group factor C\ = Cq as 
can be seen from the standard model result eq. @. The relevant Feynman graphs involving 
at least one fermion loop are those in fig. 1 with external gluons being attached at all possible 
points. We may average over the N 2 — 1 color gluons. Thus the group factors that occur are 
1/(N 2 — 1) times the group factors of the graphs obtained by adding an internal gluon line 
in all possible ways to the graphs of fig. 1, thus transforming them into four-loop graphs. 
In each of these four-loop graphs there are six gluon vertices. We may classify them by 
the resulting number of gluons attached to the outer fermion loop in fig. la-d. We ignore 
for the present quartic couplings. Each resulting graph will then have n connections on 
the outer fermion line and 6-n connections in the inner loop and we treat the graphs in 
order of decreasing n. eq. (P tells us that the group factor for a graph with any order of 
gluon attachments on the outer loop is a linear combination of the group factor from the 
graph with the gluons attached in a standard order and group factors from graphs with 
fewer gluon attachments on the outer loop. This relation may be described by a graphical 
equivalence in color space illustrated in fig. 2. Since eq. (p]) holds for all representations, the 
same equivalence is valid on a gluon or ghost loop also. Thus the linearly independent group 
factors may be found by considering planar graphs only. When gluinos are incorporated into 
the theory, each graph with a quark loop has a corresponding graph with one or more quark 
loops replaced by gluino loops. The gluino contribution to the j3 function is given by the 
quark contribution except for the replacement of the fundamental representation matrices 
by the adjoint representation matrices and by an extra factor of 1/2 for each gluino loop due 
to the Majorana nature of the gluino. This factor of 1/2 has its source in the fact that, unlike 
the case for Dirac fermions, one must ignore the direction of the gluino line in a Feynman 
graph in calculating statistical and symmetry factors. The only linearly independent group 
factor coming from graphs with n=6 may be taken as the planar graph obtained by attaching 
another gluon line at adjacent points on the outer loop of fig. lc. Using the identities of eqs. 
J|. |A2]) , the corresponding group factor is seen to be the C\ of eq. (H). 

d = — - (n f Tr (T a T a T b T b T c T c ) + ^T T (F a F a F b F b F c F c )) . (A10) 
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The traces here are easily evaluated using eqs. ([A1|) -([A"5|). Because of fig. 2, the group factor 
from the various non-planar graphs with n=6 are linear combinations of C\ and group 
factors coming from graphs with n < 6. Thus eq. ( |A10| ) can be taken to be the only linearly 
independent group factor coming from the n=6 graphs. Graphs with n=5 can obtained by 
adding a gluon line to fig. Id or fig. lc. Using fig. 2, one sees that all such graphs lead to 
group factors that are proportional to the C*2 of eq. @. 

C 2 = — (n f Tr(T a T a T b T c T d ) + -M-T T (F a F a F b F c F d )) f bcd (All) 

plus possibly a linear combination of group factors from graphs with lower n. It is clear, 
using the tracelessness of the representation matrices and the equivalence of fig. 2, that this 
is the only linearly independent group factor one can write with 5 attachments on the outer 
loop and one internal attachment. We proceed, therefore to n=4. A linearly independent 
group factor with n=4 is the C3 of eq. fllPl). 

C 3 = 4 (n f Ti(T a T b T c T d ) + r ^-Ti(F a F b F c F d ))f abe f cde . (A12) 

This appears when the extra gluon line is added connecting the gluon lines of fig. lc. A 
second group invariant corresponding to n=4 is 

(n f Tr(T a T a T b T c ) + ^Tr(F a F a F fe F c ))Tr(F b F c ). (A13) 

However, by explicit calculation this is seen to be proportional to C<i- At this point one 
should investigate the graph where to fig. Id one adds a gluon line with one leg on the outer 
loop and the other attaching to make a quartic coupling. However the group factor at a 
quartic coupling is a sum of three terms each of which is a product of two f abc . Thus in group 
space the quartic coupling is equivalent to a sum of products of triple vertices as expressed 
in fig. 3. Again this figure has the meaning that the group factor of a graph including the 
quartic coupling is a linear combination of group factors where the four legs are linked in the 
three possible ways by triple couplings. The consequence of this equivalence is that we may 
neglect all graphs with quartic couplings in determining the number and form of linearly 
independent group factors. With n=4 we also have a contribution from the two fermion 
loop topology of fig. la where a gluon line is attached at adjacent points on the outer loop 
and each of the fermion lines can represent either a quark or a gluino. The corresponding 
group factor is C4. 

C 4 = —^—(n f TiT a T a T b T c ) + r ^-Ti(F a F a F b F c )) (n f Tr(T b T c ) + ^Tr(F b F c )). (A14) 

All other graphs with four attachments on the outer loop can be seen, using fig. 2, to be 
linear combinations of C\ through C4 and possible group factors appearing at lower n. At 
n=3 we have the group factor from fig. la with the extra gluon connecting the outer and 
inner loops. The corresponding group factor is 

j^(n f Tr(T a [T b ,T<}) + ^-Tr(F a [F b , 

(n f Tr(T a T b T c ) + ^-Tr(F a F b F c )) 

= C 5 . (A15) 
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Here we have used the fact that the divergent part of fig. la when gluon crossed graphs are 
included is totally antisymmetric in a,b,c. This can be seen from the fact that the divergent 
sub-graph of fig. la with an extra gluon joining the two fermion loops is a renormalization 
of the triple gluon vertex. The other two-fermion-loop possibility at n=3 is obtained by 
connecting the extra gluon line to fig. la with one end on a triple gluon vertex. The 
corresponding group factor is proportional to C5 as can be seen by using fig. 2. The only 
remaining possibility with n=3 is 

1 - (n f Tr(T a [T b ,T c ]) + n~ g /2Tr(F a [F b ,F c })Tr(F a F b F c ), (A16) 

coming for example from adding a gluon line to fig. lb or fig. Id. However this is easily seen 
to be proportional to C3. Finally we consider graphs with n=2. One such graph coming 
from fig. la is equivalent to a graph treated earlier under interchange of the outer and inner 
loops. The only other group invariants one can construct with n=2 are of the form 

C = 1 (n f Tr(T a T b ) + ^Tr(F a F fc )) ( y A r Ti(F a F b F c F c ) + BTr(F a F e )Tr(F e F b )) (A17) 

all of which are proportional to C3. The B term, in fact, never appears in the three- loop (3 
function since after cutting a gluon line in the four-loop graph to make the external gluons, 
the graph must remain one-particle-irreducible. This exhausts the graphs contributing to 
the gluon propagator renormalization. 

The ghost propagator renormalization contributes no new group factors. This can be 
seen by noting that the ghost-ghost-gluon vertex has the same group factor, f abc , as the 
triple gluon vertex. The ghost propagator corrections can be obtained from the topologies 
of fig. 1 by adding a gluon in all possible ways which lead to at least one closed gluon loop 
and then replacing each gluon loop in turn by a ghost loop. Cutting one of the ghost lines 
then leads to a ghost propagator renormalization graph. All of these clearly have the same 
group structure as the corresponding contribution to the gluon propagator renormalization. 
In fact, since the role of the ghosts is to cancel unphysical Lorentz modes of the gluons, 
one can anticipate that the ghost propagator and vertex renormalizations contribute no new 
linearly independent group factors. It is useful, however, to see this in detail. 

We therefore consider finally contributions to the f3 function from renormalization of the 
ghost-ghost-gluon coupling. The corresponding graphs are given by attaching an external 
gluon in all possible ways to the graphs of fig. 1 and then at all possible points joined by 
a continuous gluon line attaching a pair of external ghosts and changing the continuous 
gluon line between them into a ghost line. Since the ghost-ghost-gluon vertex like the triple 
gluon vertex is proportional to f abc , the infinite corrections to the vertex must also have this 
group structure. Since f abc f abc = Ca{N 2 — 1), all the group factors from the renormalization 
of this vertex can be obtained by considering all the vacuum to vacuum, five loop graphs 
including a fermion loop and a ghost loop that remain one-particle irreducible after at 
least one particular gluon-ghost-ghost vertex is excised. The group factors that occur in 
the gluon-ghost-ghost vertex renormalization are then the group factors of these five-loop 
graphs divided by Ca{N 2 — 1). Each of these graphs consists of a ghost loop and fermion 
loop with 3 < n < 6 gluon vertices on the ghost loop. Four representative graphs are shown 
in fig. 4. The three-loop ghost-ghost-gluon vertex correction can be restored by excising 
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a ghost-ghost-gluon vertex from these five-loop graphs in all possible ways that leave the 
graph one-particle irreducible. Since there are eight gluon vertices in a five loop graph, all 
the graphs with six attachments to the ghost loop will have two attachments to the fermion 
loop as in fig. 4d. Each such graph, therefore, will have a factor 

71 ~ 71 ~ 

n f Tr(T a T b ) + -j-Ti(F a F b ) = (n f T + -j-C A )5 ab . (A18) 

The group factor for such graphs will therefore be (nfT+rigC A /2) / C A times the group factor 
(C A ) of the three-loop gluon propagator correction with no fermion loops. The resultant 
group factor is C3. A graph such as fig. 4a with three gluons attached to the ghost loop 
clearly gives no new group factors since 

-Tr(F 9 F h F i ) = —4 r. (A19) 



C A (N 2 -1) v ' 2(W 

That is, a graph with a ghost loop attached to the rest of the graph by three gluons has 
the same group factor as the graph where the ghost loop is shrunk to a point at a triple 
gluon vertex leading to one of the graphs treated in the gluon propagator renormalization. 
Graphs with 5 vertices on the ghost loop (e.g, fig. 4b) have the group factors 

\ r (n f Ti(T b T c T d ) + ^lT T (F b F c F d ))Tr(F a F a F b F c F d ) (A20) 

1 ) 2 



C A (N* - 1) 



and 



— — ^ -(rifTr(T b T c T d ) + —Ti(F b F c F d ))Ti(F a F b F a F c F d ). (A21) 

C A {N — 1) v 2 

Each of these has a group factor proportional to C3 as can be seen by using eqs. (|A1|) -( [A3| ). 
Graphs where there is a triple gluon vertex can be reduced to a linear combination of eqs. 
( |A20|) and (|A21|) using the equivalence of fig. 2. The graph of fig. 4c where the fermion and 
ghost loop are joined by four gluons has, when all the crossed gluon graphs are considered, 
the group factor 

I / 71 ~ 

/ m f^r, ^h^Qrpdx u g rp / jpa 776 t-tc jpd 



, —{n f Ti{T a T b T c T d ) + —Tr(F a F F c F ))[ATr({F a , F ]F°F ) 
C A \N — 1) ^ 2 /\ 

+BTr(F a [F b , F c ]F d ) + CTr(F a F b F c F d - F c F b F a F d )) . (A22) 



The symmetry of the divergent part of these graphs is determined from the fact that the 
divergent sub-graph is a renormalization of the quartic gluon vertex. Each of these three 
terms is proportional to C3. The group factor for fig. 4f is also proportional to C3 as is 
that for the case of fig. 4f with fermion and gluon loops interchanged. Finally there is also a 
two-fermion-loop graph renormalizing the ghost-ghost-gluon vertex. This is shown in fig. 4e 
where the vertex renormalization is obtained by excising and discarding a ghost-ghost-gluon 
vertex that leaves a connection from the ghost loop to each of the fermion loops. The group 
factor for fig. 4e is 
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Ca(n2 _ 1} (n^(n*n + ^Tr(F a F«F c )) 

(n/Tr(T c T d ) + ^Tr(F c F d ))Tr(F a F fe F d ) = (A23) 

This completes the proof that, when gluinos are included, there are only six linearly inde- 
pendent group factors in the three-loop j3 function for the SU(N) strong coupling constant, 
namely C\ and the five group factors given in eqs. (|3|)-(|l2|). This leads to the result stated 
in eqs. ([?]) and (|T3"D of the text. 
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FIG. 1. Basic topologies for the calculation of the gluon propagator renormalization. Solid 
lines represent fermions, either quarks or gluinos, and wavy lines represent gluons. The pair of 
external gluons are attached at all possible positions. 



FIG. 2. Diagrammatic equivalence in group space. The group factor from a graph containing 
the left hand side is a linear combination of the group factors from graphs in which each of the 
terms on the right hand side is inserted in its place. 



FIG. 3. Graphical equality in group space between a quartic coupling and three pairs of triple 
couplings. As in fig. 2, the equality has the meaning that the group factor of a graph containing the 
quartic coupling is a linear combination of the group factors obtained by replacing that coupling 
in turn by each of the subgraphs shown on the right hand side. 



FIG. 4. Basic topologies with a least one fermion loop for the renormalization of the 
ghost-ghost-gluon vertex. Fermions, gluons, and ghosts are indicated by solid, wavy, and dashed 
lines respectively. 
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Fig. 2 



